Coherent spin mixing dynamics in a spin-1 atomic condensate 
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We study the coherent off-equilibrium spin mixing inside an atomic condensate. Using mean field 
theory and adopting the single spatial mode approximation (SMA), the condensate spin dynamics 
is found to be well described by that of a nonrigid pendulum, and displays a variety of periodic 
oscillations in an external magnetic field. Our results illuminate several recent experimental obser- 
vations and provide critical insights into the observation of coherent interaction-driven oscillations 
in a spin-1 condensate. 
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Bose- Einstein condensation (BEC) has been one of the 
most active topics in physics for over a decade, and yet 
interest in this field remains impressively high. Recent 
experiments showcase the rich versatility of control over 
the atomic superfluid, e.g. the BEC-BCS crossover 0, 
2], quantized vortices 0, 0, 0, condensates in optical 
lattices Q, and low dimensional quantum gases 0, @. 
While most of these efforts involve condensates of atoms 
in a single Zeeman state, activities in spinor condensates 
0, have recently received significant boost with the 
addition of three new spin-1 experiments [T1IIT2I liH , Il4| . 

In a spinor condensate, atomic hyperfine spin degree 
of freedom becomes accessible with the use of a far-off 
resonant optical trap instead of a magnetic trap. For 
atoms in the F = 1 ground state manifold, the presence 
of Zeeman degeneracy and spin de pen dent atom-atom in- 
teractions |^ flfl[lll[iall^ . lr^ . ll8lll9| | leads to interesting 
condensate spin dynamics. In this article, we study spin 
mixing inside a spin-1 condensate [T7L TTflL l2Cj. focusing 
on the interaction-driven coherent oscillations within a 
mean field description. Unlike the pioneering studies on 
this subject as in Refs. 0, 0], we will highlight the 
important role of an external magnetic field, which is 
present in all experiments to date. 

Recently, a beautiful experiment has finally observed 
the long predicted Joscphson type coherent nonlinear os- 
cillations with a scalar condensate in a spatial double well 
potential pi) . 

Although spin mixing driven by the internal spin- 
dependent interaction (not of the nature of a Rabi os- 
cillation as driven by an external field (2^, HH), has 
been observed in both F = 1 and F — 2 condensates 
0- [O- El- the coherence of this process has not 
yet been investigated. Over-damped single oscillations 
in spinpopulations have been observed in earlier experi- 
ments 24] although their interpretation has been limited 
because evolution from the initial (meta-stable) states 
was noise-driven. The main experimental obstacles to 
observe more oscillations are the dissipative atomic col- 
lisions among the condensed atoms and the decoherence 
collisions with noncondensed atoms [T2I E3 | . A promis- 
ing future direction relies on increased atomic detection 
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FIG. 1: (Color online) The dependence of oscillation period 
T on magnetic field B for a 87 Rb condensate from our model 
(solid line, Eq. I8al 1. the results from a full numerical simu- 
lation without the use of SMA are denoted by (*). 



sensitivity, thus the use of smaller condensates as in the 
experiment of Ref. |2l| , with lower number densities and 
at lower temperatures, two favorable conditions for the 
single spatial mode approximation (SMA) 0, ^ . 

The initial atomic population distribution in Fig. ^ 
corresponds to the (equilibrium) ground state at a mag- 
netic field (B-field) of 0.07 Gauss and with a zero magne- 
tization (m = 0), specified by p (Q) w 0.644 and 8(0) = 
with c » 0.614 Hz (these symbols are defined later). As 
in the case of no B-ficlds [l?], EH, Il3 > the initial relative 
phases among the three components depend on the spin- 
dependent atom-atom interaction being ferromagnetic 
(0) or antiferromagnetic (tt), inside an external magnetic 
field |25j . Starting with the initial phases and population 
distributions at Bg — 0.07 Gauss, we instantaneously 
change the £?-field to a different value, the atomic con- 
densate distributions thus become off-equilibrium, and 
the coherent dynamics starts according to the mean field 
theory. Within the SMA, we find such off-equilibrium 
dynamics of a spin-1 condensate corresponds to that of 
a nonrigid pendulum, which can be characterized using 
semiclassical trajectories in the phase space. 
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Our main result is illustrated in Fig. ^ where we have 
plotted the dependence of oscillation period on the ex- 
ternal magnetic field. The parameters are close to the 
experiment |26| : where the spin- independent trap is har- 
monic V = (M /2)(uj 2 x 2 + LOy-y 2 + uo 2 z z 2 ) with u x — uj y — 
(27r)240 Hz and u> z = (27r)24 Hz. The condensate con- 
tains AT = 1, 000 87 Rb atoms with an average density of 
(n) w 1.7 x 10 13 cm" 3 . 

As illustrated in Fig. ^ the spin mixing dynamics 
within the SMA corresponds to a typical pendulum, with 
the quadratic Zeeman energy playing the important role 
of the total energy At small change of i?-field, the equiv- 
alent pendulum undergoes a small amplitude oscillation, 
approximately harmonic with a period independent of 
the energy or oscillation amplitude; increasing of the to- 
tal energy leads to a longer oscillation period as the pen- 
dulum becomes increasingly nonlinear. At a critical field 
B c , when the effective total energy is just enough to bring 
the pendulum to the completely up or top position, the 
period approaches infinity as for the homoclinic orbit of 
a pendulum; Upon further increasing the energy (or B), 
the pendulum starts to rotate around and the period be- 
comes smaller with increasing energy as the pendulum 
rotates faster and faster. 

Our system of a spin-1 atomic Bose gas inside an exter- 
nal magnetic field is described by the Hamiltonian [l5l[T^ | 



n 



dr 



4 (-|g v2 + y + E i) ^ + y^VJV^i 
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where repeated indices are summed, and ipi(r) (4) ^ s 
the field operator that annihilates (creates) an atom in 
the i-th hyperfine state (\F = l,i = +1,0,-1), here- 
after \i)) at location r. M is the mass of an atom. In- 
teraction terms with coefficients Co and ci describe re- 
spectively elastic collisions of spin-1 atoms, expressed in 
terms of the scattering length ao (a 2 ) for two spin-1 atoms 
in the combined symmetric channel of total spin (2), 
c = Airh 2 {a + 2a 2 )/3M and c 2 = ^fi 2 (a 2 - a )/3M. 
F~t=x.y,z are spin-1 matrices. Assuming the external mag- 
netic field B to be along the quantization axis (z), the 
Zeeman shift on an atom in state \i) becomes (the Breit- 
Rabi formula 27]) 



E± = 
Eo = 



— |p T giftiB - ^hfs \/l ± a 



— — Evfi?Q a/T" 
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where £hfs is hyperfine splitting, and gj is the Lande 
g-factor for an atom with nuclear spin I. pi is the nu- 
clear magneton and a = (gipiB+gjpsB)/ J?hfs with gj 
representing Lande ^-factor for a valence electron with a 
total angular momentum J. /is is the Bohr magneton. 

The field operators ipi evolve according to the Heisen- 
berg operator equation of motion. At near-zero temper- 
ature and when the total number of condensed atoms 



(N) is large, the condensate is essentially described by 
the mean field <fii — (ipi). Neglecting quantum fluctua- 
tions, they form a set of coupled Gross-Pitaevskii (GP) 
equations, from which we can simulate the mean field 
off-equilibrium dynamics more accurately at various ex- 
ternal magnetic fields without using the SMA. 

Our simplified model is based on the well-known fact 
that for both 87 Rb (ferromagnetic) and 23 Na (antiferro- 
magnetic) atoms, the spin dependent interaction oc |ca | 
is much weaker than the density dependent interaction 
oc I co I - This leads to the validity of the SMA, where 
we adopt the mode function (j>(r) as determined from 
the spin-independent part of the Hamiltonian Ti. s — 
-{Ti 2 /2M)V 2 + V + c n 0013. 

Thus we define 



N^(t)<f>(r) exp(-vrt/ft), 



(2) 



where TL s 4>{y) = pcf>(r) and J dr\<fi(r)\ 2 = 1. We arrive at 
the coupled spinor equations 



m± = e ± c± + c[{ P± + p Q - p T )e± + CoC T i 



(3) 



with c = c 2 N f dr\<j)(r) 



\S,i\ 2 - It is easy to verify 



that the total atom number and atomic magnetization 
are conserved, i.e. y\. pj = 1, p + — p- = m, and m = 
(N + — N-)/N is a constant of motion. 

We use r)=(E-- E+)/2 and 8 = (E- + E+ - 2E )/2 
to parameterize the linear and quadratic Zeeman effect. 
We further transform 

£+ -> €+exp[-i(E -r))t/h], 

Co -> ^ exp[-iE t/h], 

C_ - £-exp[-i(Eo + v)t/h] 

to eliminate the Eq and 77 dependence, and take £j = 
v /pJe _i6 ' J . After some simplification, we obtain the fol- 
lowing dynamic equations for spin mixing inside a spin-1 
condensate 



2 C r— 

Po = ft-zWl 1 



Po) 



m s'm9, 



28 2c, 

T^ 1 



2po) 

2c \ (l- j o )(l-2 i oo)-m 2 



n) 



y/(l - po) 2 - m 2 



(4) 



cos 6, (5) 



where 9 = 9+ + 9- — 26q is the relative phase. These two 
coupled equations give rise to a classical dynamics of a 
nonrigid pendulum, whose energy functional (or Hamil- 
tonian) can also be derived within the SMA as in 0] 



£ = cp 
+8(1 



(1 - Po) + pa) 2 -m 2 cos 9 

Po)- 



(6) 



It is easy to check that po 
(2/h)d£/dp . 



-(2/h)d£/d9 and 6 
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FIG. 2: Equal-energy contours for a condensate of Rb atoms 
(upper panel) with B = 0.05 Gauss, |c| = (2-7r)0.5 Hz, and 
m = 0; of 23 Na atoms (lower panel) with B — 0.015 Gauss, 
jcj = (2tt)0.5 Hz, and m = 0.3. 



The contour plot of £ in Fig. [3 displays several types 
of oscillation as in a pendulum. The dynamics of spin 
mixing described by Eqs. |@| [SJ in a magnetic field is 
conservative, as also recognized and studied numerically 
in Ref. [28|, The corresponding phase space trajectory 
is therefore confined to stay on the equal energy contour. 
Quite generally, po oscillates in a magnetic field. Rewrit- 
ing equation 0} as 



(Po) 2 = - s {[£-6(l-po)][{2cp + S){l-p )-e\ 
n 

-(cp m) 2 }, (7) 
we can compute the oscillation period according to 



T = <b —dp 
Po 



V2h g (vir~fp 

V-5c \fx 3 - Xi 



for c < 0, (8a) 



and 



T 



V2h 



X3—X2 
X3-X1 



for c > 0. 



(8b) 



K(k) is the elliptic integral of the first kind, and 
2^=1,2,3 are the roots of po = (order as x\ < 
x 2 < £3) (Fig- EJ- The period for a rigid 
pendulum, described by u + sinu = 0, is T = 
4y / 2K[ v / 2/(E + 1) }/y/E+l at an energy E > 1 and 

T = 4y/2F[arc&m(y/(E + l)/2 ), y/2/(E + 1) ]/V^TT 
when —1<E<1. Here E is the energy of the rigid 
pendulum. 

The time evolution of po can be expressed in terms of 
the Jacobian elliptic function cn(.,.), 



Po(t) = x 2 - (x 2 - xi)cn 2 70 + t\/-25c(x3 - xi), k , 
for c < 0, (9a) 
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FIG. 3: (Color online) The dependence of cubic roots Xj on 
the external magnetic field for 87 Rb atoms (left) and 23 Na 
atoms (right). Other parameters are |c| = (27r)0.5 Hz, po(0) = 
0.6, 0(0) = 0, and m = for 87 Rb; 0(0) = tt and m = 0.3 for 
23 Na. 
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FIG. 4: (Color online) The magnetic field dependence of 
the oscillation period for 87 Rb atoms (left) and 23 Na atoms 
(right). Other parameters are the same as in Fig. 



and 

Po{t) = x 3 - (23 - £ 2 )cn 2 
for c > 0, 



To 



+ ty/2Sc(x3 — Xi), k 



(9b) 



cn 2 (7 ,fc) 

,2/ 



M'2 



70 depends on the initial state, 
po(0)]/(x 2 - xi) if c < and cn 2 (7 ,/c) = [x 3 
Po(0)]/(x 3 - x 2 ) if c> 0. For 87 Rb atoms (c < 0), 70 = 
if po(0) = x\ and 70 = K(k) if p (0) = x 2 . The solutions 
of po are oscillatory between x\ and x 2 if c < (between 
x 2 and xz if c > 0), except when x 2 — X3 (x 2 = Xi 
if c > 0), where the solution becomes homoclinic, i.e., 
limt^oo po = 1 and the corresponding period is infinity 
for m — (Fig. 0}. 

We further observe from Fig. 0] that when the total 
magnetization is varied the peak of the oscillation period 
essentially stays at the same magnetic field for ferromag- 
netic interactions. The solution becomes periodic when 
m/0 since po can at most reach 1 — m. It turns out that 
the critical solution of an infinitely long oscillation period 
occurs when po(i — > 00) = 1, or equivalently £ = 0, which 
gives 8{B C ) — |c|po(l+cos 9) with p and 9 the initial con- 
ditions. At B = we reproduce the same result as in Ref. 
[T^|. The rapid decreasing of the period when B > B c 
is consistent with the recent numerical simulations by 
Schmaljohann et al. [l3|]. For antiferromagnetic interac- 
tions, however, the peak of the oscillation period shows 



a strong dependence on the magnetization, and asymp- 
totically we find po(t — * oo) = 0, i.e., £ = S(B C ) which is 
equivalent to 8(B C ) = c[(l — po) + ^/(l — po) 2 — fn 2 cos 9] . 

Substituting the solution po(t) into Eq. we can 

solve for 9(t). Furthermore we can find the time depen- 
dence of 6± and 9 through the following 



9± 



-- [d+cpo + cp < — 

n y 1 — po ±m 



cos 9] 



h 



[(1 - p ) + y/(l- po) 2 -m^cos9]. 



Finally, we consider the evolution of the averaged total 
spin. As was recently demonstrated by Higbie et al, the 
averaged spin of a condensate or its magnetization can 
be directly probed with non-destructive phase contrast 
imaging j29j. Alternatively, the magnetization dynamics 
can be inferred from component populations of a spinor 
condensate, which are directly measurable using Stern- 
Gerlach effect in an inhomogeneous magnetic field. We 
first illustrate the quadratic Zeeman effect on the spin 
dynamics of a noninteracting condensate. For a state 
£ = (£+, £o, £-) T > the total spin average is (F) = (£\F x x+ 
F y y + F z z\i) with 



(F x ) = V2Re[|£ |(|e 
(F y ) = \/2Im[M (|e 



\ e i(9o-e+) + |^_| e J («o-e_) 
|^_| e *(flo-fl-) 



+ e 



i(e -9+) 



(F, 



m. 



As an interesting case, we take the initial state as £(0) = 

iT 

V(l - Po)/2, yfpo, V(l - w/2 • po is a constant. 
We find at time t that 



(F x )+i(F y ) = 2^p (l-p )cos(5t/h)e- ir > t / n , 

(F z ) = m = 0. (10) 

It spirals toward and away from the origin in the (F x )- 
(F y ) plane. The linear Zeeman effect causes spin pro- 
cessing around the magnetic field [z axis), while the 
quadratic Zeeman effect makes spin average oscillate. 

The spin evolution becomes quite different when atom 
interaction is present. For the same initial conditions (of 
the above), the total averaged spin at time t becomes 



(F x )+i(F y ) = 2 V 'po(l-po)cos(fl/2)e- i ^/ R 5 

(F z ) = 0, (11) 

which can be conveniently confirmed from the phase 
space contour plot of Fig. [2J where 9 is confined to 
oscillate around zero for ferromagnetic interactions and 
around ir for antifcrromagnetic interactions if B < B c . 
Note that po and 9 are time-dependent for interact- 
ing condensates. Figure [3] exemplifies this oscillation 
in terms of the allowed regions (shaded) of (F x ) and 
(F y ) for interacting condensates in contrast to non- 
interacting ones. For ferromagnetic interactions, the 
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FIG. 5: Two dimensional projection of the averaged spin evo- 
lution (shaded region) for a condensate with zero magneti- 
zation of noninteracting atoms (middle), in comparison with 
atoms of ferromagnetic (left) and antiferromagnetic interac- 
tions (right). 



allowed region is defined by two radii. 



One of them, 
de- 



n = ^2 Po (0)[(i - MO)) + 7(T^o"(o)) 2 - ni 

pends on the initial condition, while the other (rg) is 
solely determined by the quadratic Zeeman effect. We 
find r B > 77 if B < B , < r B < rj if B < B < B c , 
and Tb = if B > B c . There exists a forbidden region at 
the center for a ferromagnetically interacting condensate 
if B < B c . This region shrinks to zero when B > B c . 
Exactly at B = B c , an interesting attractor-like feature 
arises and the average spin gradually spirals towards the 
origin (at the center) and becomes trapped eventually 
after an infinitely long time. For antifcrromagnetic in- 
teractions, the allowed region generally becomes smaller 
than that for a noninteracting condensate as shown in 
Fig. [SI for 77i = (or B c — 0). The radius of the shaded 
(allowed) region depends on the quadratic Zeeman effect, 
while the forbidden region approaches zero as B — > oo. 

For the general case of m 7^ 0, the al- 
lowed region is in between the two radii 

^0(0)1(1 - p (0)) ± v/(Wo(0)) 2 - m 2 ] for a 
noninteracting gas, For ferromagnetic interactions, 
the averaged spin behaves similar to the case of 
777 = considered above, except now the forbid- 
den region shrinks gradually to a minimum nonzero 

value of yj2po{0)[(l - p Q (0)) - V{l^p (0)) 2 - m 2 ] 
when B — > 00. In this case, there exists no B c 
or homoclinic orbits. For antiferromagnetic inter- 
actions, the analogous radius 773 decreases from 

ri = ^0(0)1(1 - po(0)) - ^/(Wo(0)) 2 - 777 2 ] to 
zero while B increases from zero to B c . At B — B c the 
attractor-like feature remains. When B is increased from 
B C1 tb increases from zero, and crosses 77 at B = B , 
finally approaches the radius of the allowed region for a 
non-interacting condensate when B — ► 00. 

Before concluding, we hope to make some estimates to 
support the use of the mean field theory, i.e. treating 
the atomic field operators as c-numbers. Intuitively, we 
would expect that this is a reasonable approximation as 
the total numbers of atom, at 1000, although not macro- 
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scopic, is definitely 'large'. In fact, the recent double 
well experiment that confirmed the coherent nonlinear 
Josephson oscillations of the mean field theory, is at a 
similar level of numbers of atoms |21j . A rigorous discus- 
sion of this point in terms of the quantum phase diffusions 
in a spin-1 condensate is a rather involved procedure, and 
will not be reproduced here |30l |. Instead, we illuminate 
the validity of mean field theory as follows. First, we look 
at the total atom number fluctuations. Approximating 
the spinor condensate as a one component scalar, and ne- 
glecting the internal spin mixing dynamics, its total over- 
all phase spreads after a time of r c » N/[a(N)(co(n))] 
[3l| . with cr(iV) ~ \/~N the standard deviation of the 
atom numbers from taking c-number approximations of 
the atomic field operators. In our case, this time is about 
0.2 second, short compared to a typical Josephson type 
oscillation period at ~ 1 second. We believe, however, 
this is not a critical issue as we are not studying phase 
sensitive phenomena involving the overall phase as in an 
interference experiment. Instead, we are interested here 
in the relative phase dynamics between different conden- 
sate components, whose oscillation time scale is given 
by the much smaller value of the spin-dependent inter- 
action coefficient ci\ thus we should compare the coher- 
ent classical oscillation period of ~ 1 second with the 
much longer time t' c w N/[<7(N)(c 2 (n))} [T^. ~ 50 sec- 
onds (for 87 Rb). This then leads to a favorable condition 
for adopting the mean field theory in our study. Alter- 



natively, we can reach the same conclusion from a direct 
investigation of the oscillation period T Eq. I|8al) , which 
contains a simple N dependence cx l/\/jV. We find that 
\T(N ± y/N) - T(N)\/T{N) = l/(2y/N), is only about 
2%, indicating the overall validity of the mean field the- 
ory. 

In conclusion we have studied the off-equilibrium in- 
teraction driven collective oscillations inside an atomic 
condensate in an external uniform magnetic field. The 
dynamics of spin mixing is found to be well described 
by a nonrigid pendulum due to the conservation of atom 
numbers and atomic magnetization. In particular, we 
find that there exists an interesting class of critical tra- 
jectories whose oscillation periods approach infinity. Our 
study illuminates the use of quadratic Zeeman shift to 
probe pendulum-like oscillations in a spin-1 condensate 
and provides the complete spin mixing dynamics ana- 
lytically. It provides the much needed theoretical guid- 
ance for the eventual experimental detection of coherent 
macroscopic oscillations in a spinor condensate. 

Note added in proof: We have recently observed many 
of the coherent oscillatory behavior discussed in this work 
and have submitted a publication describing these exper- 
iments. 

W. Zhang thanks Dr. Su Yi for several discussions of 
this project. This work is supported by NSF and NASA. 
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